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1. Introduction 

The infrared properties of the Landau-gauge QCD Green functions have been triggering many 
studies in the last few years, mainly involving both lattice (see for instance Refs. [1-3] ) and 
continuum approaches (see for instance Refs. [4-12] using Dyson-Schwinger equations (DSE), 
[13-16] using the refined Gribov-Zwanziger formalism, [17] using the Curci-Ferrari model as an 
effective description) or based on the infrared mapping of A0 4 and Yang-Mills teories [18]. 

Most of the DSE analysis take advantage of approximating the ghost-gluon vertex by a con- 
stant when truncating the infinite tower of the relevant equations. To be more precise, only the 
behaviour of the involved transverse form factor needs to be approximated by a constant for the 
purpose of truncating the ghost propagator DSE (see also the recent paper [19]). Most of the ammo 
for the approximation of the ghost-gluon vertex and GPDSE truncation is mainly provided by the 
Taylor theorem which is widely known as a non-renormalization one [20]. This claims that, in 
the particular kinematical configurations defined by a vanishing incoming ghost-momentum, no 
non-zero radiative correction survives for the Landau-gauge ghost-gluon vertex, which takes thus 
its tree-level expression at any perturbative order [20]. This statement entails two important conse- 
quences: that (i) the bare ghost-gluon vertex, defined by 



is UV-finite for any kinematical configuration; and that (ii) in the specific MOM scheme where 
the renormalization point is taken with a vanishing incoming ghost momentum (which we named 
the Taylor scheme -T-scheme- in [21]) the renormalization constant for the ghost-gluon vertex 1 is 
exactly 1. This allows for a precise lattice determination of the strong running coupling defined 
in T-scheme 2 that can be confronted to perturbative prediction in order to estimate [21, 23, 
24]. The transverse character of the gluon propagator, in the Landau gauge GPDSE kernel, allows 
to project out the fy-iorm factor from expression ( |1.1[ ). Therefore, we have chosen to call the 
surviving piece, i.e. H\, the transverse form factor. Extending this nomenclature, we call H2 the 
longitudinal vertex form factor. 

A precise determination of the ghost dressing function obtained by solving the GPDSE, to 
be confronted for instance with lattice data, thus requires a correspondingly precise knowledge 
of the transverse form factor H\ . To investigate the non-perturbative structure of that form factor 
was one of main goals of ref. [25], which we will pay attention to in this contribution. That was 
accomplished, as will be also shown here, by studying the non-perturbative OPE corrections to the 
perturbative form factor. The non- vanishing dimension-two gluon condensate, (A 2 ), plays a crucial 
role for the leading non-perturbative contribution. This condensate was proved not to be neglegible 
when studying the running of the QCD Green functions [21, 23, 26-31], but very specially from the 
analysis of the running of the T-scheme coupling computed by means of lattice simulations with 
2 + 1 + 1 twisted-mass dynamical flavours in ref. [32], as will be also shown here. 

1 It can also be straightforwardly concluded that the Landau-gauge ghost-gluon vertex renormalization constant, Zi , 
is exactly 1 in the MS scheme. 

2 It has been also used to define an effective charge that could be applied for phenomenological purposes [22]. 
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2. The OPE for the T-scheme running coupling and (A 2 ) 

A very recent analysis [32] of the strong coupling in T-scheme, 

= 1 {Taylor theor. } 

a T {q 2 ) = Urn G(q 2 ,A 2 )F 2 (q 2 ,A 2 )(H l (q,0)+H 2 (q,0)) , (2.1) 

A->0 47T 

where the bare ghost and gluon dressing functions, G and F, were computed from lattice simula- 
tions with 2 + 1 + 1 twisted-mass dynamical flavours, provided with a very strong evidence for the 
necessity of applying non-perturbative corrections to account for the running of the coupling, 



a T {ix z ) = a£ ert (ju^ 



V 



1+^ K >W(^) 4( ^_ 1} 



(2.2) 



where R(a,GCo), y§~ and (Xj n can be obtained in perturbation and through the OPE analysis (see 
ref. [32]). This can be in Fig. [j], borrowed from [32]. The left plot shows the lattice data for 
the Taylor coupling multiplied by the square of the momentum plotted in terms of the four-loop 
perturbative value of the coupling at the same momentum, with a best-fit of Aj^g- which gives 

«ms(" j z) = 0.1198(9) , (2.3) 

in very good agreement with PDG [33]. One should notice that the departure from zero for the 
lattice data in the plot can be only explained by non-perturbative contributions and that the Wilson 
coefficient for the Landau-gauge gluon condensate successfully describes the non-flat behaviour 



from the data up to a pert ~ 0.32 (the solid line given by Eq. (2.2)), which roughly corresponds to 



p ~ 3.5 GeV. The right plot shows the best-fit of Eq. ( |2.2| ) to the lattice data, with a best-fit value 
of 

8 2 (ql)(A 2 ) R , ql = 4.5(4) GeV 2 , (2.4) 

where = 100 GeV 2 . These are of course estimates for A^g- and the gluon condensate including 
2 degenerate light quarks and two non-degenerates heavier ones. This allows to run in perturbation 
the coupling value up to cross safely the bottom mass threshold and run again up to get a "physical" 
estimate for the coupling at the Z°-mass scale, that can be directly (and successfully) compared 
with experimental results. The same analysis had been previously applied to lattice data in pure 
Yang-Mills QCD (Nf = 0) [21] and with Nf = 2 twisted-mass dynamical flavours [23], providing 
with estimates for the gluon condensate slightly depending on the number of light quark flavours, 
although roughly ranging from 4-6 GeV 2 , when statistical uncertainties are included. 
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Figure 1: (Left) deviation from the lattice data with respect to the prediction of the four-loop perturbative 
theory, with the best-fit of A7-, plotted in terms of the perturbative running; the solid line shows the leading 
non-perturbative OPE prediction, Eq. (2.2)'s r.h.s. (Right) The strong running coupling in Taylor scheme 
defined by Eq. (2.1) obtained over a large momentum range from lattice QCD in ref. [32]. The dotted 
line stands here for the best-fit with Eq. (2.2), while the solid one includes a higher-order power correction 
effectively behaving as ~ 1//A 

3. The OPE for the ghost-gluon vertex 

We will now pay attention to (and describe with more details) the OPE procedure applied in 
ref. [25] for the Landau-gauge ghost-gluon vertex, 

V abc { _ qk . q _ k) = Yf>J { _ qk . q _ k) G bb^ q _ k)F aa' {q)F c^ k) 

= J d 4 y d 4 x e i{q - k) - x e iky ( T (c c (y)A* (jc^O)) ^ (3.1) 
is quite similar to the one described in ref. [21]. Here, the OPE expansion shall read 



V^(-q,k;q-k) = (d )f c (q,k) 

+ id 2 )ffi/{q,k) (:A*;,(0)At',(0):) + 
where do accounts for the purely perturbative contribution to the vertex, while 

= 2 /[ 1 ] + 2/iV2/l 2 ]+4/ [3l +2/W+/ [51 



where 3 



w 1 



q-k *\ 




Nc 2 {q-k)aqa bc 
2 8 q 2 (q-kf Vtre ^ [ q,k ' q ] ' 



Nc is the number of colours. 



(3.2) 



(3.3) 



(3.4) 
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the tree-level ghost-gluon vertex being 



while 



V t f e U-q,k;q-k) 



k 2 q 2 {q-k) 2 




q-k 



-V-2--- £i. q - 



/W q 
I k- 



(3.5) 



Nc 2 (q-k)qk a bc 
2 8 k 2 (q-k) 2 q,k,q h 



and 



C 



Nc i^qayabc ( _ ata _ k) 

4 * ^2^2 v tvss,n\ q^q K ) ■ 



The other contributions, /P 4 5 !, can be immediately written using the previous results: 



Z [3] 



q- k fo 




V 

Ncf Nc8 2 
4 k 2 4 



y / 



and 




2/W+/M = 2: 



I + 




(3-6) 



(3.7) 



(3.8) 



(3.9) 



In all cases, the blue bubble refers to a contraction of the colour and Lorentz indices with \8 a bgp.v 



We have also introduced the notation g^ v (^) = g^ v — -*jr for the transverse projector. 
Then, one obtains 



w 



f = g 2 I s v (q,k)+s F (k)+s F (q)+s G (k-q) j V^-q^q-k) 



(3.10) 
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with 



sc(q) = s P {q) =^y , 
T 

s V {q,K) 2 ^^ 2 ^_^ )2 , ~ k 2( q _ k y > k 2 q 2 



2 (k-q)-k k q 



(3.11) 



if (resp. sg) comes from the OPE corrections to the external ghost (resp. gluon) propagator, 



i.e. from the non-proper diagrams in Eqs. (3.8,3-9) and sy from the proper vertex correction in 



Eqs. (3.4,3.7). The OPE corrections to the proper vertex can be obtained from there, remembering 
that the bare ghost-gluon vertex reads as: 

rf(-q,k;q-k) = -gof abc (q^Hi(q,k) + (q - k)^H 2 (q,k)^j . (3.12) 

Then, in Landau gauge, only the form factor H\ survives in the Green function to give: 

V^ hc (-q,k;q-k) = -ig Q f abc q^g^^q-k)Hi{q,k) G((q-k) 2 ) F(q 2 ) F{k 2 ) , (3.13) 

where G and F are the gluon and ghost dressing functions for which we previously computed the 
non-perturbative OPE corrections. Thus, one would have: 



Hi{q,k) = Hf* t {q,k)[ l+s v {q,k) 



(A 2 ) 



4(A£-1 



+ 



4 q-\k-\q- 2 k- 2 ) 



(3.14) 



All the non-proper corrections have been removed in the usual way. Now, in order to compare 
all over the available range of momenta with some current lattice data, we need to try to model 
the low-momentum behaviour, non-accessible via our previous OPE analysis. This is also done in 
ref. [25] and will be described in the next section. 



4. The Euclidean ghost-gluon vertex and the present lattice data 

We will now devote this section to model the transverse form factor of the Euclidean ghost- 



gluon vertex, Hi, on the ground provided by Eq. ( |3.14| ), and compare the result with some recent 
lattice data for this form factor computed in different kinematical configurations. 

4.1 The model for the Euclidean ghost-gluon vertex 

In Euclidean metrics the bare ghost-gluon vertex can be written very similarly as: 



(4.1) 



K^A-q^,q-k) = ig f bc [q^Hiiq^ + iq-k)^^) 



where the form factor H\ plays a crucial role when solving the ghost-propagator Dyson-Schwinger 
equation (GPDSE) as discussed above. The OPE non-perturbative corrections to the form factor 



H\ were obtained in Eq. (3.14), but that result is in principle only reliable for large enough k,q 
and q — k since the SVZ factorization on which it relies might not be valid for low momenta. 



Nevertheless, in the following we will propose a very simple conjecture to extend Eq. (3. 14) to any 
momenta, which will provide us with a calculational model for the ghost-gluon vertex to continue 



6 



The dimension-two gluon condensate 



J. Rodriguez-Quintero 



research with. In particular, the perturbative part of the ghost-gluon vertex is usually approximated 
by a constant behaviour 4 ; if we then apply a finite renormalization prescription such that 



Z x (jl z )H x {q,k)\. 



1 



(4.2) 



where the renormalization momentum, /I 2 , for a given kinematical configuration (for instance, 
q — k = and q 2 = k 2 = /x 2 ) is chosen to be large enough, on the basis of Eq. ( |3.14| ), we can 
conjecture that 



1 + 



N c g z (A 2 ) 
8(A£-1) 



\/ k 2 q 2 cos 6 
x [ -^-^ — + 2 



q 2 — \J k 2 q 2 cos 6 



k 2 q 2 + m 



IR 



q 2 (q 2 +k 2 - lyfq 2 ! 2 cos d) + 



m 



IR 



k 2 — \/k 2 q 2 cos 6 



k 2 {q 2 + k 2 - l^fq 2 !? cos d) + 



in 



IR/ 



(4.3) 



gives a reasonable description of the ghost-gluon form factor Hi all over the range of its momenta 



q and k, where 6 stands for the angle between them. The purpose of Eq. (4.3) is to keep the main 
features of the momentum behaviour of the ghost-gluon form factor provided by the OPE analysis 
and to give a rough description of the deep IR only with the introduction of some IR mass scale, 
miR, that can be thought as some sort of IR regulator (related to some kind of effective gluon mass) 
and which is mainly aimed to avoid the spurious singularities resulting from the OPE expansion in 
momentum inverse powers. 

Of course, \i 2 being large enough, Z\ can be approximated by some constant value for any \i 2 
because the logarithmic behaviour of the ghost-gluon in perturbation theory has been proven to be 



very smooth [36, 37]. Thus Eq. (4.3) provides us with a very economical model because one needs 
nothing but some infrared mass parameter which, being related to the gluon mass, is expected to 
be ~ 1 GeV, to parametrize the deep infrared behaviour of the ghost-gluon transverse form factor. 
It should be noted that, although the ghost-gluon vertex form factors do not diverge, the gluon 
condensate does: it is the product of the Wilson coefficient and of the condensate which is expected 
to remain finite thanks to a delicate compensation of singularities, in the sum-rules approach. Then, 
both the Wilson coefficient and the condensate should be renormalized, by fixing some particular 
prescription, and both shall depend on a chosen renormalization momentum. This is explained in 
detail in the work of ref. [31] where the Wilson coefficient for the quark propagator is included 
at order &(a 4 ). Since we carry out the computation of the Wilson coefficient in Eqs. (3.4-3.9) 



only at tree-level and neglect its logarithmic dependence, we do not specify any renormalization 
momentum for the gluon condensate. 



At least for large momenta, this seems to be the case in lattice simulations [34, 35] for several kinematical config- 
urations, and so is confirmed also by the perturbative calculations in Refs [36, 37]. 
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4.2 Comparing with available lattice data 



Some SU(2) and SU(3) lattice results for the ghost-gluon vertex are available in the literature 
(see for instance Refs. [34, 35]). In particular, Cucchieri et al. have computed the tree-level- 
tensor form factor in SU(2) for three kinematics (p = q — k is the gluon momentum and (p the 
angle between the gluon and ghost momenta) : (1) p 2 = q 2 and (p = 7t/2, (2) p 2 = and (3) 
<p = jt/3 with p 2 = q 2 = k 2 ([34]). As can be seen in the left pane of Fig. ||, when we plug 
g 2 {A 2 } = 6 GeV 2 , miR = 1.4 GeV and Zf 1 = 1.04 into Eq. (fD|), the prediction for the form factor 
Hi appears to agree pretty well with the SU(2) lattice data for the three kinematical configurations 
above mentioned. It should be also noted that the OPE prediction without the infrared completion 



through the introduction of wiir given by Eqs. ( p.llU3.14j ), plotted with dotted lines, also accounts 
very well for lattice data in the intermediate momenta region, above roughly 2.5 GeV. Furthermore, 
the value for the gluon condensate lies in the same ballpark as the estimates 5 obtained from the 
SU(3) analysis of the Taylor coupling in ref. [21]. This strongly supports that the OPE analysis 
indeed captures the kinematical structure for the form factor H\ . On the other hand, the infrared 
mass scale, as supposed, is of the order of 1 GeV and the (very close to 1) value for Z^ 1 accounts 
reasonably for perturbative value of H\ that we approximated by a constant. Thus, after paying 
the economical price of incorporating only these two more parameters, both taking also reasonable 
values, we are left with a reliable closed formula for the form factor at any gluon or ghost momenta. 
This is a very useful ingredient for the numerical integration of the GPDSE. 



X '■' 
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— qW; 


felt/3 


— k 2 =0 






5 
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Figure 2: (left) The predictions of the model for H\ (solid lines) in Eq. (4.3) and the OPE result without 
infrared completion (dotted lines) in Eq. (3.14) confronted to the SU(2) lattice data borrowed from ref. [34] 
for the three kinematical configurations described in the main text: (k — q) 2 = q 2 for (p = n/3 (violet) and 
(p = n/2 (red) and (k — q) 2 — (green), (right) The results of the model for Hi in the SU(3) case and for the 
following kinematical configurations: k 2 = q 2 for 9 = n/3 (violet) and 9 — n/2 (red) and k 2 = (green). 

The SU(3) results published by the authors of ref. [35] appear to be very noisy and cannot 



be invoked to properly discriminate whether a constant behaviour close to 1 or Eq. (4.3) accounts 
better for them. However, we can now take the mass parameter, mrR, from the previous SU(2) 
analysis and the well-known SU(3) value for g 2 (A 2 ), assume Z^ 1 = 1 and use Eq. (43) to predict 



5 In ref. [21], g 2 (A 2 ) is evaluated through an OPE formula with a tree-level Wilson coefficient to be 5.1 GeV 2 at a 
renormalization point of 10 GeV. 
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the ghost-gluon transverse form factor, H\. This is shown in the right plot of in Fig. ^ for three 
different kinematical configurations and, as can be seen, the deviations from 1 appear to be very 
small in all the cases and compatible with the results shown in Fig. 4 of ref. [35] for the vanishing 
gluon momentum case. 

4.3 The Taylor kinematics and the asymmetric gluon-ghost vertex 

The kinematical configuarition where the incoming ghost momentum goes to zero (k — > 0) is 
especially interesting, as this is the one which one defines the T-scheme for (as shown in sec.||). 
Let us pay some additional attention to it. In this particular kinematical limit, it is exactly obtained 
up to all perturbative 6 orders: 

rfe M (-9,0tf) = -gf bc q, • (4.4) 



At any order of perturbation theory, this implies that//] (q,0) +H2(q,0) = 1. According to Eq. (4.3), 
one would have 

HlM )=z r '( 1+ iV CJ gil J -^ s ). (4.5) 

An interesting question to investigate is whether such a genuine non-perturbative correction still 
survives for the full ghost-gluon vertex in the Taylor limit, and not only for the transverse form 
factor H{. Said otherwise, does eq.([0]) maintains its validity upon adding the OPE corrections 
related to (A 2 ) to it, or not? This question have been properly addressed in ref. [25], where the 
asymmetric ghost-gluon vertex is studied and one is left with 

rf c (-q,e;q-e) = -gf abc [q^ H x {q,e) + (q-e)^H 2 (q,e)^ +■■■ (4.6) 



where 



H 1 (q,e)=Hr(q,e)+N c g Aq - £) - q A " 



q 2 {q-e) 2 4(JVg-l) 

2 { {q-e)-qf (A 2 

q 2 {q-ef 4(A£-1) ' 



H 2 (q,e) = Hr(q,e) ~ N c g 2 "«~ c ^j 77 ^ TT . (4.7) 



Thus, after taking the limit e —¥ 0, one will have: 



Ff c (-q,0;q) = sf bc q, ( ^M) + H^O) \ = -gf bc q, . (4.8) 



Thus, the result applied in Eq. (2.2) to define the T-scheme coupling is recovered. This confirms 
that no non-perturbative OPE correction survives in the proper ghost-gluon vertex in the Taylor 
limit, although both form factors Hi and H 2 separately undergo such a kind of correction. 



6 The same argument of Taylor's perturbative proof still works if one consider the Landau-gauge ghost-gluon vertex 
DSE: a vanishing ghost momentum entering in the vertex implies the contraction of the gluon-momentum transversal 
projector with the gluon momentum itself for any dressed diagram. Thus, Taylor's theorem is still in order within the 
non-perturbative DSE framework. Some attention have been recently paid to the ghost-gluon vertex DSE [38]. 
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5. Conclusions 

The momentum behaviour of the ghost-gluon vertex has been studied in the framework of the 
operator product expansion (OPE). This approach has already been proved in the literature to be 
very fruitful in describing the running of two-point gluon and quark Green functions and of the 
strong coupling computed in several MOM renormalization schemes. In particular, the coupling 
derived from the ghost-gluon vertex renormalized in T-scheme, which can be directly computed 
from nothing else but the ghost and gluon propagators, led to a very accurate determination of 
Aqcd when its running with momenta obtained from the lattice estimate has been confronted with 
the OPE prediction, although only after accounting properly for a dimension-two gluon condensate, 
(A ). The same approach is followed to study the ghost-gluon vertex and special attention is payed 
to the transverse form factor for the ghost-gluon vertex, precisely the one that plays a crucial role 
for the truncation and resolution of the GPDSE and that, in this context, is usually approximated by 
a constant. In this way, a genuine non-perturbative OPE collection to the perturbative part for the 
transverse form factor is obtained. A very simple conjecture, made to extend the OPE description 
beyond the momentum range where the SVZ factorization is supposed to work, left us with a 
simple model describing the momentum behaviour of the ghost-gluon form factor. This model was 
shown to agree pretty well with a lattice SU(2) computation for several kinematical configurations 
of the ghost-gluon vertex, when the value for the gluon condensate is in the same ballpark as the 
one that is estimated from the running of the T-scheme coupling in SU(3) lattice gauge theory. 
This successful comparison with the available lattice data provides us with strong indications that 
(i) the OPE framework helps to account for the kinematical structure of the ghost-gluon vertex at 
intermediate momentum domain and (ii) that it helps to cook up a reliable model, providing a good 
parameterization also for the IR domain, to be plugged into the GPDSE to reproduce the ghost 
propagator lattice data. Finally, we also proved that, in the particular kinematical limit associated 
with the Taylor theorem and hence by the T-scheme (a vanishing incoming ghost momentum), the 
corrections to the longitudinal and transverse form factors cancel against each other and that the 
tree-level result for the ghost-gluon vertex is thus recovered. 

Acknowledgements: J. R-Q acknowledges the Spanish MICINN for the support by the re- 
search project FPA2009-10773 and "Junta de Andalucia" by P07FQM02962. 
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